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Abstract 

Bethe-Salpeter equation is applied to nucleon-nucleon elastic scat¬ 
tering at the intermediate energy. The differential cross section and 
the polarization are calculated in terms of the phase shift analysis 
method using the two-body potential derived from the Bethe-Salpeter 
equation. The lowest-order Born approximation for the K-matrix is 
corrected by including the inverse square part of the potential. 

One of the main purpose for study of nuclei is to understand the nuclear 
force between nucleons and the effects on the nuclear many-body system. The 
meson exchange picture and the two-body interaction play a decisive role to 
describe the nuclear structure. To give an example the higher-order dia¬ 
grams beyond the Hartree or the Hartree-Fock approximation elucidate why 
the lowest-order treatments of the proper self-energy reproduce observables 
of the nuclear system|l]. And the nuclear effective interactions as illustrated 
by the pairing interaction in the nonrelativistic formulation are largely con¬ 
nected with the two-body correlations in nuclei. 

Strictly the meson exchange potential used to calculate founds on the 
screened coulomb type, thus, the relation with two-body interactions in free 
space is still an open question. In order to treat the relativistic system com¬ 
posed of two nucleons the Bethe-Salpeter (BS) equation is indispensable to 
the underlying framework of the theory and the validity on the held theoret¬ 
ical points of view has been reported in the literature for a long period. The 
geometry of the space-time is four dimensional and the interaction contains 
relativistic effects such as the retardation effect, which is often neglected to 
make the calculation tractable in a lot of models of the potential. 

The integral equation on the two-particle Green function is used to derive 
the BS equation. By equating the pole terms on the both sides the differen¬ 
tial equation is derived for the bound state composed of two nucleons, that 
is, deuteron. The singular potential is necessary to evaluate the quantity and 
the BS equation is expected to provide analyses at the short-range or the high 
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momentum transfer region. Recently it has been shown that the binding en¬ 
ergy and the electric quadrupole moment of deuteron are obtained by solving 
the ladder approximated BS equation with the asymptotic approximation!^. 
It is the main subject of the present study to investigate scattering states 
between two nucleons, for which the boundary condition exists in the asymp¬ 
totic region and determines the observables. 

The BS amplitude is expanded by a set of the Gamma matrices Tj 
and which results in the simultaneous equations on the expansion coeffi¬ 
cients. They act as the wave functions Tj {i=S,V,PS,AV,T), in which each 
subscript denotes the scalar, vector, pseudoscalar, axial-vector and anti¬ 
symmetric tensor sectors respectively. By the Fourier transform they are 
converted from the four dimensional momentum space to the conhguration 
space ones so as to proceed calculations for the present study. 

When terms of the total energy are taken into account exactly the equation 
of the vector sector Ty is independent of the other components provided that 
the auxiliary condition Po^o(p) = 0 is imposed on the zeroth vector compo¬ 
nent To(p), where po is the zeroth component of the relative four-momentum 
in the center of mass system. Calculating the binding energy of deuteron 
the vector component is removed because the zeroth component is assigned 
the spin S'=0. The three polar components of the anti-symmetric tensor are 
substituted to construct the S'=l bound state of deuteron. 

Doing the Gamma matrix expansion and applying the auxiliary condition 
to the origin of the relative time (t=0) in the conhguration space the equation 
results in the form equivalent to the Schrodinger eigenvalue equation as 
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where /i and E denote the reduced mass of two nucleons and the laboratory 
energy of the incident nucleon respectively. The dummy index is applicable 
to j. The subscripts i,j=V, AV or PS on 5=0 and i,j=x, y, z the cartesian 
coordinate of the polar-tensor on 5=1. When 5=0 the potential is assumed 
to be Vji{x) = 6jiVi{x). It has been reported on the study of the electric 
quadrupole moment of deuteron that the mixing term in the polar-tensor 
equation nearly cancelled out the joint effects of the axial-vector and the 
pseudoscalar components, thus, it is neglected in the following calculation 
for 5=1 part as well as that of the pseudoscalar equation applied for 5=0 in 
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the numerical calculations. 

It is remarkable that Eq. (1) is derived without the nonrelativistic ap¬ 
proximation as the expansion by n/c, then, may take any positive value 
of the incident energy for scattering states. We make a suggestion that the 
pseudoscalar component \hp 5 (a:;) is used for the singlet (*S'=0) state to con¬ 
struct the proton-neutron [p-n) elastic scattering along with a set of 'by (*) 
{i = x,y,z) for the triplet (5=1) state which has also been applied to de¬ 
scribe the bound state quantities. Since the mixing term in the pseudoscalar 
equation is dropped out these two spin states 5=0 and 5=1 do not couple 
mutually similar to the two-body interaction potential Vi(x). 

When one investigates the scattering problem the partial wave analysis 
method is generally used for proton-proton (p-p) case and formulated by the 
M-matrix with the isospin T=T3=lj3]. This formulation is extended to p-n 
elastic scattering, in which both of the isosinglet (T=0) and the isotriplet 
(T=l) matrices are prepared to obtain the resulting observables. There 
is reason why the zeroth vector component is not adapted for 5=0, that is, 
the spin function Ty is symmetric ^Ty = Ty, then, the parity (—1)^ = -l-(—)1 
are assigned T = 0(1) states accordingly. On the other hand, the spin 5=1 
is carried by Tp of the tensor components which is symmetric ^Tp = Tp and 
the assignment of the parity is equal to that in the usual formulation. 

For the mixing between 5=0 and 5=1 is not assumed as stated above, 
the partial waves are decomposed into three parts, (5,L, J) = (0,L,L), (1,L,L) 
and (1 ,J=f 1,J). The last part allows mixing between two states L=J=f 1 and 
which is represented by three real parameters 5jpij and ej. The nuclear bar 
phase shift method is used to express them|3]. The mixing parameters ej 
play a decisive role in giving suitable magnitudes on the M-matrix elements 
Mio and Mqi by which the value of the polarization may be influenced much. 

In order to determine the phase shift parameters for p-n elastic scattering 
the formulas 

POO 

tan = —2A:p / r'^dr jL(kr)Vi(r)'ijjL(kr), (3) 

Jo 

are applied with use of the potential derived from the BS equation in 5=0 
case. Here, jiihr) and 'ipLikr) are the spherical Bessel function of the order 
L and the exact solution under the potential in Eq. (1) respectively. The 
5=1 case is analogous to Eq. (3). It corresponds to the Born approximation 
in the K-matrix formulation when jL(fcr) replaces 'ipiikr) by the lowest-order 
approximation on the strength of the potential Vi{r). The K-matrix is 2x2 
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real symmetric matrix and related to the phase shift parameters of the S- 
matrix represented by the eqnation S' = (1 + — iK) in the algebraic 

notation. 

In the present meson exchange model the isoscalar scalar meson a, the 
isoscalar vector meson u, the isovector psendoscalar meson vr and the isovec¬ 
tor vector meson p are taken into acconnt. Valnes of the conpling constants 
and masses are determined appropriately referring to calcnlations for the rel¬ 
ativistic nnclear many-body system and the strnctnre of denteron[2]. The 
electromagnetic interaction is tnrned off in the present stndy which has an 
inflnence on the resnlts particnlarly at the forward direction of the center of 
mass system {Of. < 10°) in p-p system. 

Concerning the interaction of pion the psendovector conpling is employed 
and the propagator is mnltiplied by the fonr dimensional cnt-off fnnction 
AV(A2 - p2) in the momentnm space. By means of the cnt-off procednre 
the inverse fonrth power potential shape (~ r“^) at r —)■ 0 is modihed to the 
inverse sqnare potential one (~ r“^) and which makes calculations for the 
phase shifts (Eq. (3)) feasible. Another important feature of the pseudovec¬ 
tor coupling of the pion interaction is the property of the angular momentum 
changing (L —)■ L ± 2) in the two-body potential which effectively acts like 
the tensor force giving a suitable size of the mixing parameter unlike the case 
of the psendoscalar coupling interaction. 

While the Born approximation works well for the Coulomb potential, the 
nuclear potential is much stronger at the short-range region and for the exact 
wave function 'ipiihr) we need to take into account the inverse square po¬ 
tential originating in the Feynman propagator therein. Instead of the exact 
solution by numerical procedure only the leading-order inverse square part of 
the potential is left and the approximate solution ipL{kr, g) is substituted for 
'ipL{kr) in the present study. For scattering states Eq. (1) under the inverse 
square potential is solved analytically and it is given as 


'^hi.kr^g) 


sec(- j- 
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The order of the Bessel function u = [L + — g is dependent on the 

strength g of the inverse square potential part {V{r) ~ -g M~^ where 
M is the nucleon mass. When we solve the equation the Neumann function 
part is dropped by assuming the phase shift 5l = (L-1-1/2 — z/)7r/2 on the 
inverse square potential. Thus, the elements of the K-matrix is improved by 
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the approximate function '4’L{kr,g) in this manner. It is verihed that the 
spin S'=0 part is not influenced much by the procedure so the correction is 
not done for the Born term. Concerning the spin S=1 part we correct only 
for (L, J) = (0,1) state as to the p-n scattering case to enhance the phase shift 
parameters of the J=1 part in the K-matrix suitably. 

Instead of the direct expansion of the K-matrix Eq. (3) is used for the 
(0,1) element by multiplying the factor F{A) to the Born term as 


<L'\V\iy> = <L'\V\L> F{A), (5) 

F{A) = <L'\ V{A) \ u> / <L' \ V{A) \L> . (6) 

Here, the bracket notation is used to clarify some points on the cut-off A for 
the propagator of pion explained previously. The argument A denotes that 
the cut-off A in the potential is allowed to vary as a parameter from A^ to 
the inhnity oo. On the other hand, the hxed cut-off A^ is in the potential 
V and u to make the integral of the Born term be convergent. The value of 
To ~ 500 - 600 MeV is tentatively used by taking account of the calculation 
for the binding energy of deuteron[2]. As a result of the numerical calculation 
it is shown that F{A) remaines hnite as T—)-cxd and the results with F{oo) 
give remarkable effects on the scattering phenomena compared with F{Ao). 
So, F{oo) is chosen to calculate the higher-order corrections in the present 
study. 

The differential cross section for p-n elastic scattering is given as 

= - y TrM^M^t^ (7) 

dQ 8 j’To 1 

where and denote the M-matrices of the isosinglet and isotriplet re¬ 
spectively. In Fig. 1 the angular dependence of the differential cross section 
at the laboratory energy E = 310 MeV is shown in two cases (V)the lowest- 
order approximation in the K-matrix theory and (K)inclusion of effects of 
the higher-order terms by the procedure stated above. By including the 
corrections the curve becomes deeper at the center of mass scattering angle 
6c = 90° and shows the desirable V-shape. Thus, the short-range character 
of the nuclear force is responsible for the angular dependence largely and 
in fact contributes to determine the incident energy dependence of the p-n 
scattering favorably. 

The polarization is another means of investigating properties of the nu¬ 
clear force and the accurate phase shift parameters with the differential cross 
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section equally. The results of the calculation for polarization are shown in 
Fig. 2 where the signs (V) and (K) denote same as those in Fig. 1. Since the 
present formulation assumes the identical two fermions in the isospin space 
the curve is symmetric at 6c = 90° different from the asymmetry seen in the 
experiments. The result of the lowest-order calculation (V) is not satisfying 
for lack of the sufficient tensor force. Then, by including the higher-order 
effects the curve turns the sign successfully. The higher-order calculation 
results in an underestimate particularly at the forward direction (6*c < 90°) 
compared with the experiment in which the maximum value reaches 0.4 at 
around 9c = 30° |1]. The discrepancy may be attributed to the approximate 
form of 'ipL{kr,g) adopted for calculating the K-matrix, which would be im¬ 
proved by performing calculations next to leading order perturbatively or 
correcting the other (T, J) states in addition to (0,1) state. 

The BS formalism is also applicable to p-p elastic scattering by adding the 
inhomogeneous term to the BS equation j5j. While the role of the various com¬ 
ponents in the BS amplitude is not necessarily conhrmed, the pseudoscalar 
component is assigned the 5=0 scattering state as well as the p-n case. The 
S-wave (T=0) is essential for p-p elastic scattering to interpret the observed 
straight line shape at 6c > 10° of the differential cross section. As to the 
5=1 state a set of the polar-tensor components is appropriate following the 
p-n case. 

The differential cross section for p-p elastic scattering is given by chang¬ 
ing the sum over T in Eq. (7) as 2 X]t=i- If is verihed that 

the numerical calculation results in an overestimate of the experimental one 
when parameters of the coupling constants, the meson masses and the cut¬ 
off parameter of pion in the present model are same as those of the p-n 
elastic scattering. In order to make up for the ladder approximation the 
pion-nucleon pseudovector coupling constant / is reduced to / ~0.05 from 
the standard value / ~ 1[6] at E = 310 MeV incident energy. The suppression 
of the neutral vr® meson exchange force indicates that the lowest-order ladder 
diagram in the irreducible kernel is not sufficient and prompts us to correct 
by the higher-order multi pion processes largely. 

The triplet P states (L,J) = (1,0), (1,1) and (1,2) are fundamental to con¬ 
struct p-p elastic scattering and among them only (1,2) state is related to the 
mixing parameter 62. Then, the (1,2) element of the 2x2 K-matrix is calcu¬ 
lated by using 'ipi{kr, g') instead of i\{kr)^ in which g' denotes the strength of 
the inverse square part of the potential. One complicated issue in the present 
formulation is that the pion exchange force is attractive when T=1 in the 
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S=1 state, therefore, the g' becomes too strong to leave the z/' real even in 
L=l. The imaginary part invalidates the phase shift method to describe the 
p-p elastic scattering consequently. Dealing with the situation we make a 
decision to move the parameter g' at —oo < g' <\,hY which the results are 
made feasible to compare with the experiments maintaining the framework 
of the phase shift method. 

It has been found that the required strong repulsive force is in S'=T=0 of 
the p-n system. The inverse square potential in the pseudoscalar equation is 
mainly given by the isoscalar a and u mesons as 

^ (27r)2 {27i)‘^ml. 

Here, ga^ and g^^ denote the coupling constants accordingly. Meanwhile, the 
leading order of the pion exchange force with the cut-off function is attrac¬ 
tive and less important. Ultimately the expected strong repulsive force arises 
from the u meson exchange interaction. 

The replacement for g is interpreted as manifestation of ir^-p-n three- 
body system in the intermediate state. Since the positive pion (tt^) carries 
the isospin T=1 the intermediate p-n state is recognized as T=0 provided 
that the transition from the initial p-p state is done by a charge-independent 
interaction conserving the isospin. By the anti-symmetric property of the 
pseudoscalar sector (^T ps = —T ps) giving the spin function on S'=0 the or¬ 
bital angular momentum of the p-n system comes to L=1 and followed by 
P-wave (T=l) state of to construct the odd parity J=2 state as a whole. 

The differential cross section for p-p elastic scattering is shown in Fig. 3. 
While the calculation (V) is done by the lowest-order approximation the pion- 
nucleon pseudovector coupling constant / is changed as the above mentioned 
way to improve on the ladder approximated irreducible kernel effectively. As 
seen in the calculation (K) the weakened tt® exchange force is compensated 
by the higher-order correction which improves the shape of the curve well 
in comparison with the calculation (V) keeping the strong angular depen¬ 
dence under the Born approximation. The procedure is effective against the 
underestimate of the polarization that is larger than the present calculation 
seen in Fig. 4[4,7]. Particularly, enhancement of the mixing parameter €2 by 
the exact treatment of the K-matrix is probably significant to reproduce the 
polarization like the p-n elastic scattering case. 

In the p-p system it has been found that the coupling constant of pion 
decreases as the incident energy of proton increases in order to reproduce the 
differential cross section of the elastic scattering. The energy dependence of 


7 



the p-p system is more obvious than the p-n system in which there is little 
need to adjust the parameter between the 10 MeV and 200 MeV. The BS for¬ 
malism is suitable to comprehend two-body nucleon-nucleon system for once 
the interacting lagrangian is set up the higher-order calculations are carried 
out systematically about each component of the BS amplitude by various 
methods in the quantum mechanics. Determining which one is appropriate 
to describe each channel of the scattering phenomena experimental facts on 
the spin observables are useful. 
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Figure Captions 


Figure 1: Proton-neutron differential cross section as a function of the center 
of mass scattering angle at the laboratory energy of 310 MeV. 

V:The result of the calculation by the Born approximation. 

Kdncluding the higher-order corrections. 

Figure 2: Proton-neutron polarization as a function of the center of mass scat¬ 
tering angle at the laboratory energy of 310 MeV. The signs V and K denote same 
as those in Fig. 1. The experimental data is from ref. [1] for comparison. 

Figure 3: Proton-proton differential cross section as a function of the center of 
mass scattering angle at the laboratory energy of 310 MeV. 

V:The result of the calculation by the Born approximation. 

Kdncluding the higher-order corrections. 

Figure 4: Proton-proton polarization as a function of the center of mass scat¬ 
tering angle at the laboratory energy of 310 MeV. The signs V and K denote same 
as those in Fig. 3. The experimental data are from ref. [1] and from ref. [7] at 
312 MeV for comparison. 
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